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1. ABSTRACT 


rectl^J f ? ( T l Cal ° r l03d) buckHng were established for orthotropic 

thTlTl P r e t , f °T edgG COnditions b * usin 6 the Rayleigh-Ritz method of 

3Xr5imen? n ? £ ^ ° f & StrUCtUral SyStem ’ ^dimensional buckling interaction curves 

sandwich na 1 ^ mteractlon surfaces were constructed for high-temperature honeycomlvcore 

andw,ch pane s supported under four different edge conditions. The interaction surfaces provide overall 

Streng ! hS and tHe d ° mainS ° f S * mmetricaI and antisymmetrical buckling 
sociated with the different edge conditions. In addition, thermal buckling curves of these sandwich 

pane s are presented. The thermal buckling conditions for the cases with and without thermal moments 
were found to be identical for the small deformation theory. 

!t” d ,” i t, Pan , elS ’ J the el l“ l ° f lransveree shear is large, and by neglecting the transverse shear 

Sing !t™«i, overpredicted considerably. Clamping the edges could greatly increase 

buckling strength, more in compression than in shear. 


2. NOMENCLATURE 


Amr i> A kt 

Fourier coefficients of trial function for w , in. 

A\j 

extensional stiffnesses of sandwich panel, lb/in. 

a 

length of sandwich panel, in. 


Qo 

edge length of square sandwich panel, in. 


a mnkl 

coefficients of characteristic equations 



Fourier coefficients of trial function for j xz , 

in/in. 

b 

width of sandwich panel, in. 


Cm ri) Cfct 

Fourier coefficients of trial function for j yZ} 

in/in. 

Dij 

bending stiffnesses of sandwich panel, in-lb 


D Qx } Dq v 

. transverse shear stiffnesses in xz,yz planes, Ib/in. 

D* 

flexural stiffness parameters, y/TJuThi, in-lb 

dx, dy 

differentials of x and y, in. 


E X , Ey 

Young’s moduli of face sheets, lb/in 2 


Ecx y Ec yi Ecx 

effective Young’s moduli of honeycomb core. 

, lb/in : 


Fmn Fourier coefficients for Mj, (in-lb)/in 

Gcxy, G C xz , G Cyz effective shear moduli of honeycomb core, lb/in 2 
^xy shear modulus of face sheets, lb/in 2 

Hrnn Fourier coefficients for My, (in-lb)/in. 

h depth of sandwich panel = distance between middle plane of two face 

sheets, in. 


depth of honeycomb core, h c = h-t s , in. 

moment of inertia, per unit width, of a face sheet taken with respect to 
horizontal centroidal axis of the sandwich panel, 

L = \tsh 2 + in 4 /in. 



I 

3 

k 

kx > 


ro X y 

kx ) ky 


ro X y 

l 

Mx) My 
M X y 

AfJ, Mj , Mjy 
m 

N x ,N y 

N X y 

^I^ N y> N Iy 

n 

Qxi Qy 
Smn 

T 

T a 

Ter 

t. 

V 

Vi 

v 2 

AVi 

AV 2 

AV 

u,v,w 

x,y,z 

a x> a y > &xy 

7xz>7yz 

c 


index, 1, 2, 3, • • • 
index, 1,2,3, ■ • - 
index, 1, 2, 3, • • • 

i N x a 2 u — NyQ? 

compressive buckling load factors in x- and y-directions, k x > v tc^D* ’ 


for a — constant 


NxyOt 


.2 _ 


constant 


shear buckling factor, k xy = for o = constant 

modified compressive buckling load factors in x- and direct ions, 

k x = = k x |, k y = ^$ = k vi for ab = a o = constant 

modified shear buckling load factor, k xy — Tp* = f° r a b 

index, 1, 2, 3, • 

bending moment intensities, (in-lb)/in. 

twisting moment intensity, (in-lb) /in. 

thermal moments, (in-lb)/in. 

number of buckle half waves in x-direction 

normal stress resultants, lb/ in. 

shear stress resultant, lb/in. 

thermal forces, lb/ in. 

number of buckle half waves in y-direction 
transverse shear force intensities, lb/ in. 

Fourier coefficients for (in-lb)/in. 

temperature, °F 

assumed temperature, °F 

critical buckling temperature, °F 

thickness of sandwich face sheets, in. 

total potential energy of sandwich panel, in-lb 

strain energy of sandwich panel, in-lb 

work done by external forces, in-lb 

component of Vi associated with a particular indicial condition, in-lb 
component of V 2 associated with a particular indicial condition, in-lb 
component of V, AV = AVi + AV 2 , in-lb 
displacement components in x-, y-, and z-direction, in. 
rectangular Cartesian coordinates 
coefficients of thermal expansion, in/in-°F 
transverse shear strains in xz- and yz-plane, in/in. 
numerical coefficient of Ny in a 

numerical factor in buckling equation, which changes with the edge condition 
numerical coefficient of iVj in a^ nfc j 
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v X y, v yz Poisson ratios of face sheets, also used for those of sandwich panel 

vCxy, vcyz) v Cxz Poisson ratios of honeycomb core 

Ptx specific weight of titanium material, lb/in 3 

Phc specific weight of titanium honeycomb core, lb/in 3 

3. INTRODUCTION 


Structural components of hypersonic flight vehicles (e.g., spacecraft, rockets, reentry vehicles, hyper- 
sonic aircraft, etc.) are subjected to hyper-thermal loadings due to hostile aerodynamic heating during as- 
cent and reentry, or due to solar radiation during spaceflights. The structural components of those vehicles 
ave to operate at elevated temperatures and are, therefore, called hot structures. Because of nonuniform 
heating (which is magnified by the cooler substructural frames which act as heat sinks) and the mechanical 
s ructural constraints, severe thermal stresses could build up in those hot structures. Excess thermal load- 
mg may induce (1) material degradation, (2) thermal creep, (3) thermal yielding, (4) thermal buckling, 
5) therma crack fracture after cool-down, etc. Any disruption of surface smoothness of the structures 
(e.g., metallic thermal protection system (ref. 1), hypersonic aircraft engine inlet structures (refs. 2 3) etc ) 
caused by the previously mentioned failure modes, especially thermal buckling, could disturb the flow’ field, 
creating hot spots which could cause serious consequences on the structures. Thus, the thermal load does 
play a key factor in the design of the hot structures. Reference 1 discusses various design concepts of hot 
and cryogenic structural components for the hypersonic flight vehicles. The potential candidates of high- 
uc ing-strength-hot-structural panels (fabricated with super alloys) for hypersonic aircraft applications 
are tubular panels, beaded panels, truss-core sandwich panels, hat-stiffened panels, honeycomt>core sand- 

TrLT V tC ' ( w' 4 ’ 5) ' The combined - loa d buckling behavior of the tubular panels was extensively 
studied by K ° et al. (ref. 4), theoretically and experimentally. The compressive buckling characteristics of 
the beaded panels were investigated by Siegel (ref. 5). 

Recently Ko and Jackson (ref. 6), and Percy and Fields (ref. 7) studied the compressive buckling 
e avior of the hat-stiffened panel designed for application to a hypersonic aircraft fuselage skin panel 
furthermore, Ko and Jackson conducted simple analysis of thermal behavior (thermal buckling of a face 
sheet) of the honeycom^core sandwich panel (ref. 8), and compared the relative combined-load buckling 
strengths of the truss-core and honeycomt>core sandwich panel (ref. 9). They also investigated the effect 
of fiber orientation of the metal-matrix face sheet on the combined-load buckling strength of honeycomb 
core sandwich panels (refs. 10, 11). Most of the past mechanical buckling analyses of the sandwich panels 
(re s. , refs. 9 12) and flat plates (refs. 13, 14) were conducted for simply supported edge conditions 
because the analysis was mathematically less involved. For the case of clamped edge conditions, one can 
ate the work by Green and Hearmon (ref. 15), who studied combined loading stability of plywood plates, 
nd Smith (ref 16), who considered only pure shear buckling of the plywood plates. Kuenzi, Erickson, 
and Zahn (ref 17) considered also shear stability of flat panels of sandwich construction. The works cited 
ere ignored the transverse shear effect in their analyses. King (ref. 18) analyzed the stability of clamped 
rectangular sandwich plates subjected to in-plane combined loadings, taking into account the rotational 
e ect of the sandwich core. He used a less-compact displacement function (which could be reduced to 
a simp er Green and Hearmon displacement function (ref. 15)), resulting in a complicated expression for 
the potential energy of the sandwich system. Most of the past thermal buckling analysis was done on 
single plates (refs. 19-22) or laminated composite plates (refs. 23-27), for which the transverse shear effect 
may e neglected. In actual application of the hot structural panels, most panel boundary conditions 
are closer to the clamped edges rather than to the simply supported edges. Therefore, this report will 
consider the combined-load mechanical and thermal buckling of sandwich panels under different types of 
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edge conditions by taking into account the transverse shear effect. The report also compares the buckling 
interaction curves and surfaces for different edge conditions. 

4. DESCRIPTION OF THE PROBLEM 

Figure 1 shows the geometry of a rectangular honeycomb-core sandwich panel having identical face 
sheets. The extensional and bending stiffnesses of the sandwich panel will be provided by the two face 
sheets only, and the transverse shear stiffnesses by the honeycomb core only. 

This type of sandwich panel, when fabricated with a high-temperature alloy (e.g., titanium), becomes 
the so-called hot structure, and could be a potential candidate for hypersonic aircraft structural applications 
(ref. 1). Figure 2 shows the sandwich panel subjected to combined compressive and shear loadings in its 
middle plane. The conventional Rayleigh-Ritz method of minimizing the panel’s total potential energy will 
be used in the combined-load buckling analysis, accounting for the transverse shear effect (fig. 3). The 
sandwich panel will be supported under four different edge conditions (fig. 4) 

• Case 1: Four edges simply supported (4S edge condition), 

• Case 2: Four edges clamped (4C edge condition), 

• Case 3: Two sides clamped, two ends simply supported (2C2S edge condition), and 

• Case 4: Two sides simply supported, two ends clamped (2S2C edge condition) 

where sides and ends are parallel to x and y axes respectively. 

The problem is to study the effects of the panel edge condition and the panel aspect ratio on the 
combined-load buckling behavior of the sandwich panel. Case 1 has already been solved and has been 
published in reference 9. However, for completeness, some key equations for Case 1 will be repeated in this 
report. 

5. GOVERNING EQUATIONS 

5.1 Constitutive Equations 

Following the classic orthotropic thick plate theory, which accounts for the transverse shear effect, 
the membrane force intensities {N x ,N y ,N xy }, the transverse shear force intensities {Q x ,Q y }, and the 
moment intensities { M x ,M y , M xy ) in an orthotropic sandwich panel may be related to the muidle^surface 
displacement components {u, v, w ] , the transverse shear strains {7z2>7yz}> thermal forces {N x , N y , N xy ], 
and thermal moments {M x , , Mj y } through the following constitutive equations (fig. 5) 


N x 



A\i 

0 


' du 
dx 



N? 1 


= 

A 21 

A-22 

0 


du 

dy 


~ 


N xy 


0 

0 

-466 


du , 

di + 

dv 
dx . 


- "Iy ■ 
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( 2 ) 





d ( dw _ \ 



' M x ' 


£*ii £*12 0 


dx Jxz ) 


■ Ml ' 

My 

= 

£*2 1 £*22 o 


d f dw \ 

^ J 

— 

K 

M X y 


o 0 Dee 




M t 




dy l c?x ^ xz ) dx l dy ^yzj 


- 1 - 1 


’ Q x ' 


Dq x 

0 

7xz 

Qy 


0 

D Qy 

lyz 


For the sandwich panel whose extensional and bending stiffnesses are provided only by the two identical 
face sheets, and the transverse shear stiffnesses only by the honeycomb core, the extensional and the 
bending stiffness {Aij,Dij} in equations (1) and (2), and the transverse shear stiffnesses {Dq x , Dq v } in 
equation (3) may be written as 


E x 

f — VxyVyx 


Mi , £*n 


Vyx E x 

^12 , D \2 


1 ““ VxyVyx 

M\ , £*21 

= [2 t a , 2 I s \ 

VxyEy 

-422 > £*22 


1 ~ VxyVyx 

-466 j £*66 




^ ~yVyx 

G xy 


Dq x 


= [he] 


Gcxz 

Gcyz 


(4) 


(5) 


In equations (4) and (5), {E x , E y , G xy , v xy , u yx } are the elastic constants of the face sheets, 

{G c XZl Gcyz} are the effective transverse shear moduli of the sandwich core, t s is the face sheet thickness, 

h c is the sandwich core depth, and I s is the moment of inertia of each face sheet taken with respect to 

horizontal centroidal axis, given by j _ 1 f l 2 , 1 .3 

1 s — t- j*2 l s (bj 


where h is the depth of the sandwich panel (fig. 1). The 2 in front of {£ 5 , I s } in equation (4) is associated 
with two identical face sheets. 

The thermal forces {N^ , Ny , Nj y } and the thermal moments { M Mj } appearing in equations 
(1) and (2) are defined by 


■ AfJ 

M x ' 

2 


M y 

= £ 

- N?y , 

< M xy - 

t=i 

L 


[ t s T i} (-i yU^T t 


1 v xyV yx 

t VxyVyx 

0 


VyxEj r 




xyVyx 


-1 VxyVy : r 

0 


0 

Gxy 


Qfj I 

a v 1 

a xy ■ 


(7) 


where i — 1 , 2 are respectively associated with the lower and the upper face sheets, {a x , a y , ot X y} 
are the coefficients of thermal expansion of the face sheet material, and [ ]» (i = 1, 2) implies that the 
material properties are associated with temperature 7* (i = 1,2). The thermal force and thermal moment 
contributions from the honeycomb core were neglected. 
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5.2 Energy Equations 


Based on the small deformation theory, the strain energy V\ of the heated sandwich panel may be 
written as (refs. 23, 24, 26, 27) 


An / 9u \ ^ 
2 

+ J 4l2 

(?) (?) + 
\dx) \dy ) 

A 22 
2 

to 

+ 

^66 

2 

/ du 
\dy 

au\ 2 

+ ^J 


9 

_dx 

/ dw 

\dF ~ 

‘ 7xz) 

] 2 + £>12 

d 

dx 

( dw 

■)][: 

L( 

dw N 

"a 7yz 

D22 

' d 

/ dw 

\ 

1 2 £*66 

d 

( dw 

\ . 

d_ { 

c au> 

2 

.dy 

\d^~ 

’7»a J 

J + — 

.dy 


dx * 

^ ■ 7 " 

, D Qx 
2 

7 1 , + 

7 2 - 

lyz 

<(?) 

mrT ( du 
1 ” Nxy \dy + 

di) 

_ n t f<b) 
y lay/ 


' d 

( dw 

\ 

1 . .T r 

d ( 

dw > 




+M X 

dx 

\dx ~ 

Ixz J 


9y ~ 7v *y 

•j 




' d 

( dw 

\ 

d ( dw 

Ml . 




+ MFy 

dy 

(a® “ 

-7xzJ 

dx V dy 

J j j dxdy 




( 8 ) 


For the buckling problem, the work done V 2 by the in-plane forces to produce transverse deflection is 
given by 

2 /Q. A /Q*.A / \ 21 


The total potential energy V of the sandwich panel is then 

y = vi + V 2 


dxdy 


( 9 ) 


( 10 ) 


For pure mechanical buckling problems, the strain energy equation (8) reduces to 

r b (Dn\dfdw VI 2 , _ r d (dw \]\d (dw_^ \ 

+ ° 12 [di (to " lxz )\ [dy lay 7yz )\ 


Vi 


-Ff 


+ 


D 22 


■ a 

f dw \ 

l 2 , A>6 

a 

( dw \ 

. a 

( dw V 

ay 


J + 2 

.dy 

\dF~ lxz ) 

+ 9i 

7y ‘) J 


T 2 


+ £«£ 7 2 + 

' O ixz ' 


£0y 2 \ 
2 lyz J 


dxdy 


(ii) 


and the signs of the in- plane forces {N x , N y} N xy } in equation (9) are to be reversed because the loadings 
are in the negative direction according to the sign convention shown in fig. 5. Namely, 


N t 


— AT*, A/x * Afyi A^x y 


-N, 


xy 


( 12 ) 


The sign for the shear force N xy is immaterial because, in the plate buckling problems in which shear 
is present, the eigenvalues (or shear buckling loads) always occur in pairs which are equal in magnitude 
but opposite in sign. 

For thermal buckling problems, if the sandwich panel is under uniform temperature for which 
{A/J, , M^y} = 0, ( symmetrical temperature distribution in the panel depth direction will also produce 
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zero and « the panel edges are restrained against lateral and in-plane displacements 

"r {N - w are uniform **“■ the “ 

N * = -"Z. N, = -Nl j\r„ = -A£ (13) 

proWem is *— 10 the ^ 

{M 't 1 mT diftre " ce betwe “ ‘ he *»» ft* sheets, then the thermal moments 

iU pmblem Aslm he s T ^ ‘ h f Pr ° b ' em iS "° lonser “ eigenvalue problem but a bend- 

the conditions for nnh tmn & T’ ° U ‘i r & ° S ° VC *° r U '' "' It ‘ "tv*) ' n terms of the thermal moments, and 
conditions for unbounded values of {», 7 „, ly ,) wi n give , ho buckling |oads 

5.3 Panel Boundary Conditions 

edge?co^itas h P "* 1 “ *° ^ SUPP ° rted “ four “*«“ ‘ b ' ^-g f<™ cases of boundary (or 

For mechanical buckling: 

Case 1. Four edges simply supported (4S edge condition) 

x = 0,o: w = M x = lyz = 0 ( 14) 

y = 0,6: w = M y = lxz = 0 ( 15) 

Case 2. Four edges clamped (4C edge condition) 

x = 0,a : w = || = 7xz = lyz = 0 (i 6 ) 

y = 0,6: w = |^ = 7xz = 7yz = 0 (17) 

Case 3. Two sides clamped, two ends simply supported (2C2S edge condition) 

x = 0,a: w = M x = 7yz = 0 ( lg j 

y = 0,6: w = || = lxz = 7j/2 = 0 ( 19 ) 

Case 4. Two sides simply supported, two ends clamped (2S2C edge condition) 

x = 0,o : w = = lxz = lyz = o ( 20 ) 

y = °' b: w = M y = lxz = 0 (21) 

For thermal buckling: 

In addition to the above boundary conditions, the following edge condition is to be imposed 

3 ? 0, (2 ; ti — v ~ 0 (22) 

V 0, b . a ~ v — 0 /oq\ 
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6. BUCKLING ANALYSIS 


The conventional Raleigh- Ritz method of minimization of total potential energy (V) will be used in 
the buckling analysis. To use this method, one has to assume deformation functions for the sandwich 
panel in infinite series forms containing unknown coefficients. By minimizing V with respect to each of 
those unknown coefficients, one will obtain a set of simultaneous homogeneous characteristic equations for 
eigenvalue solutions (buckling loads). 

6.1 Panel Deformation Functions 

For an eigenvalue solution via the Rayleigh-Ritz method, the trial functions for the sandwich panel 
deformation {u>, 7xz, 7yz}> satisfying the boundary conditions (eqs. (14) through (21)), may be expressed 
in the following double Fourier series for different edge conditions. 

Case 1. Four edges simply supported (4S edge condition) (ref. 9) 

(24) 

(25) 

(26) 


(27) 


(28) 


(29) 

Case 3. Two sides clamped, two ends simply supported (2C2S edge condition) (ref. 15) 


w(x,y) = 

OO OO 

= sin X X ^ mn sin 

mnx 

a 

sin 

nny 

b 

(30) 


° m=ln=l 





II 

£ 

OO OO 

= sin ^ X X ^ mn cos 

mnx 

a 

sin 

nny 

b 

(31) 


® m=ln=l 




lyz(x,y ) = 

OO OO 

= cos IT X X c ' mn sin 

mnx 

a 

sin 

nny 

b 



b m=ln=l 





OO OO 

+ sin X s ^ nc>mn 

. mnx 
sm 

rt 

nny 
cos — — 
b 

(32) 


0 m=ln=l 








w 


, . r— » r— » . . Tnnx . nny 

(x,y) = X X — 8m — 

m=l n— 1 


. \ v- V- n rn ? XX . nny 

lxz(x t y)= X X 5 mn cos — sin — 

m=l n — I 

/ \ ^ ~ . mnx nny 

lyz (z, y) = X X Cmn sm — COS — 

m— 1 n — 1 


Case 2. Four edges clamped (4C edge condition) (ref. 15) 

7 rx m , . mnx . nn 

w(x, y) - sm — sm — X X Amn sin ~~ Y sm ~~b 

a m=ln=l 


/ x 71 x • 71 y r> 

7 xz(x, y) = cos — sin — 2^ X Bmn sin 

fl w 

m 

7 xx . ny 


77i7rx . nny 
— sin 


m = 1 n=l 
oo oo 


a b 

mnx , nny 


-I- sin — sin , / / mB mn cos sin 

a b i a b 

m — 1 n — I 


'ly- 


. nx Try v- ■ m7rx ■ ni[ y 

.(x, y) = sm — cosj 2L C m „ sin sm — r~ 

“ b m=ln=l 


a 


nx . ny ^ ^ „ . mnx nny 

+ sin — sm — / X nCmn sm cos 

a 6 a b 
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Case 4. Two sides simply supported, two ends clamped (2S2C edge condition) (ref. 15) 


OO OO 


W 


(x, y)= Sin — ^ ^ A r — - WX n7ry 


„ / .-‘mriSin sin — 

° m=ln=l O b 


( 33 ) 


OO OO 


7xz(x,y) = cos ~ £ B ">n sin — — sin 

“ U 1 ..I “ h 

7 TT 00 00 

+ Sin T- £ E m5 mnCOS „„ 

a a b 


Trinx , Tl'KV 
sin 


•W*,v) = «„S 53c„ nSin !2^c«2^ 

Q- „ n A 


m=l n— 1 

OO OO 


m~ 1 7i — l 


(34) 

(35) 


assum^^riaulnc^in^ 01181 ^ 1 ^ mr1 ’ &nd Cmn ^ the undetermined Fourier coefficients of the 

“JTJZ1Z ** ^ rCSPeCtiVely ’ ^ ^ » are the buckIe ^ wave numbers in 

6.2 Uniform Temperature (Zero Thermal Moments) 

6.2.1 Mechanical and Thermal Bucklings 

thermal diitorZ” Tto Th “ "1 eXten “ 1 l<>adS ' bUt ‘ S the of strained 

th :r 

If the olr rT er ' m , mcch , am l cal b “ ckli "*. ‘ he »‘omal loads are held constant during bucklihg 

thermal buckling analysis. C ° nV<in '° na methods ° f structural stability analysis could be applied to the 

ing londSdTsrtW fa' firS ‘ H 1 tHe m f “ i “!. bUCkli " g Under the C ° mbined load - 
plied directly to the therm w f ' 7^ tmg mechanical buckling equation could then be ap- 
temped “ ‘^ ‘^l^ingol the sandwich panels with constraint edges under uniform pand 

u = o = L o) *' " N *> = “<»} <«>- (13)), (MJ. «T, M T) = 0, 

6.2.2 Rayleigh-Ritz Method 

It, (T(r !mi t f11 de ^r‘ i0n r UnCtiOM ^ ' (24> ‘ hrOUgh (35 » into ener ^ equations for 
, : ", , V2 ^ q * i 51 ^ 5 of forcing functions reversed according to ea (12) or Q ft 

tilts ar: prl^“„l B eK C ° nd,t, ° nS **"“ >— «“»*-. The 

with'^p^t^to^each Ftm^r^cwfficiwit 8 ^ ^ ^ ^ V 

mn, B mnt and C mn according to the Rayleigh-Ritz principle 

dV dV dV 

dCZ: = ° (36) 


dA r 


dB. 


mn 
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there results three homogeneous simultaneous equations (i.e., characteristic equations) for each mdicial set 


of {m, n} 


oo °o f r d * /7r\ 2 

EE -iU+i^U) 

fc t v. L 


Akt + a mnkfBlcf + Qmnfc*Cfc< j 


OO OO OQ 1 

HYl [°mnl W A kt + O^nkt^ki + a mnkt C kt\ 

k l 


= 0 

= 0 


OO OO 1 

EE u + ■&«<?«] = o 

fc i 


(37) 

(38) 

(39) 


where the coefficients a% nkt (i,j = 1, 2, 3) are defined in Appendix C for different edge conditions under 
particular indicial conditions, (a£ nkt are nonzero only when {k,£} are related to given {m,n} as shown 
in Appendix C), {k x ,k y } (contained in the coefficient a“ nfc/ , Appendix C) and k xy are respectively e 
compressive and shear buckling load factors defined as 


fcjc 


N x a 2 _ N y a 2 _ NxyQ 2 
TD?' v ~ *D* ’ xy t r£>* 


(40) 


and the flexural stiffness parameter D* is defined as 

D * = ^E Li s 

1 Vxy^yx 


(41) 


In equations (37) through (39), the indices {k,£} for nonzero a% nkt terms are determined respectively 
from the given indices {m, n} as shown in Appendix C, and in equation (37), 77 is a numerical parameter 
and 6 mnki is a special delta function which is nonzero only under the indicial conditions m±k - odd and 
n±£ = odd. Both rj and 6 mn kl respectively change their numerical values and functional forms with the 

change of panel edge condition as follows. 

Case 1. 4S edge condition: 


V 

6mnki 


= 32 

= „ gwfcL — m ; m ± k = odd, n ± £ = odd 

(m 2 - fc 2 )(n 2 - F) ’ 


(42) 


Case 2. 4C edge condition: 

_ ( 16) 3 

mrtfclfm 2 + fc 2 — 2j[n 2 + 1 2 — 2] . (43) 

tirnnkt = _ ^)[( m + 0 - 4Jl(m - k) 2 - 4][(n + 1) 7 - 4)l(n - £) - 4] ’ 

m ± k = odd, n ± £ = odd 


Case 3. 2C2S edge condition: 

77 = 8 3 

c mnk£[ 2 — (n. +£ )] ; m ± fc = odd, n ± £ = odd 

6mnki ~ (m 2 - FXn 7 - £“)[(" + - 4][(n - *)* - 4] 

Case 4. 2S2C edge condition: 

7 7 = 8 3 

I o _ mnfc£[2 - (m 2 + fc 2 )] • m ± fc = odd, n ± £ = odd 

[ ° mnfc * “ (m 2 - fc 2 )(7a 2 - Fj[(m + A:)'' 4 - 4]f(m - A:) 2 - 4] 


(44) 


(45) 
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6.2.3 Characteristic Equations in Terms of Load Factors 

FVom equations (38) and (39), for each set of [mnkl\ Ih, and r >, , J 

t nriKi), ts kt and C ke may be solved in terms of A ke 


as 


B kt = 

a mnke a mnkt ~ a mnki a mnkt ^ 
C ke = ^nk( a mnk( ~ a mntf«l'n|.f . 


(46) 

(47) 




EE 


k l l ** xy 


M, 


mnkt 


4- (5i 


mnkt 


A k t = 0 


(48) 


for every integral value of {m, n}. The stiffness- 
denned as 


geometry parameter M 


mn« appearing in equation (48) is 


Mmnkt = Tjfp-(n) 


a 11 

~mnkt 
classical 
L thin plate 
theory term 


(49) 


+ 




/*2 2 /y33 93 T9 

_mnkt a mnkt ~ nke a mnkt 

transverse shear effect terms 

combined^' tekH„ 8 g '^IT the l “ d f-*”- (*-. *», W, and is suited for 

rewntten in terms of temperature as shown in the next s,a,™ 8 eq “ ati ° n (48 > needs b « 

equ^lons^s^^ed'^vjti^t^e whofe range of indicia! °c aa ,' n 8 V'^ e uu^ber of simultaneous homogeneous 
for each set of {ra , „) value5 ,. ° f md ' clal «™b, natrons of {m, n) (one infinite series equation 

- 4 u^o^ 

and {n = f orn-f = 2) because of it, ine nrst ^rm requires the indicia] restrictions {m = k or m - k = 21 

C). Thus, if m±n is even, then k±i is alsTevln and ^ ( ^ C ' 49) ’ A PP endix 

term of equation (48), the nonzero 6 mnk( requires the indTctl “ 7, ’ ^ iS &ls ° ° dd ' In the **and 

follows that (m±k) ± ( n ±i) = ( m ±n) ±7fc±^ L ^ conditions: m±* = odd andn±f = odd. It 

the same characteristics as the first term N^ely if rn±n i US ’ ° f equation ( 48 ) also has 

is odd, then k ± t is also odd. Therefore there if no rot 1 T”’ ^ k± * ™ &ls ° even ’ Likewise, if m±n 

and the odd case (antisymmetric buckling) in the si mutt 17 ^ ^ CaSC ( symmetric buckling) 

equation (48). Thus, those simultaneous homogenel u Tiua7o7 f Seated from 

are independent of each other; one group in whfch m ± 7 7 d ' V,ded int ° two which 

odd (refs. 9, 13, 14). * P ^ m ± n 15 even ’ and other group in which m ± n is 
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For the deflection coefficients A M to have nontrivial solutions under the ass lg ned values of {k x k y a }, 
the determinant of coefficients of unknown A kt of the simultaneous homogeneous equations generated rom 
equation (48) must vanish. The largest eigenvalue 1 fk xy thus obtained will give the lowest shear buckling 
load factor k xy for given {k x , k v ,^}. When the transverse shear effect is neglected (eq. (48)), {fc x , k y , k xy } 
are a function only of J, and independent of panel size. However, they will become panel-size-dependent 
if the transverse shear effect is considered. 

In the eigenvalue computations, the infinite number of the simultaneous equations and the infinite series 
summation of each equation may be truncated up to a certain identical finite number if the convergency 
of the eigenvalue solutions has reached a specified criterion for convergency. 

The determinants in terms of the coefficients of the simultaneous equations written out from equation 
(48) up to order 12 are given in Appendix D for the cases m ± n = even (symmetric buckling) 
m ± n = odd (antisymmetric buckling) for different edge conditions The determinants of older 12 were 
found to give sufficiently accurate eigenvalue solutions (ref. 9). In Appendix D, one notices that for the 
4S edge condition only, the term (eq. (48)) forms the diagonal terms of the determinant, and the 

nonzero off-diagonal terms consist of known numerical values. However, for other edge conditions, 

not only forms the diagonal terms of the determinant, but also appears in the off-diagonal terms (mixed 

with known numerical terms). 

6.2.4 Characteristic Equations in Terms of Temperature 

For thermal buckling, the main objective is to find the buckling temperature, T cr . Therefore, equation 
(48) needs to be rewritten in terms of temperature rather than load factors. For the uniform temperature 
case, the thermal forces have the following forms 



2 

1 = (Ana x + A\2Qt y ) T 

(50) 

"» T = k * D ‘ (z. 

t2 

J = (^21 <%x + ^22 <*y) T 

(51) 

N *V = k xyD* (~ 

— AQftQtxyT 

(52) 


which were obtained from equations (4) and (7) setting T\ T 2 T . 

The coefficient a“„ u appearing in equation (49) contains thermal forcing terms (Appendix C). Thus, 
a mnJc£ ma Y be written in two parts as 

(53) 


a“ nW = T&ua + [e(m, k)Nj + C(n,*K] 


where a 11 kf is the first part of o“ nfe< without the thermal forcing terms, £(m, k) and C(«» 0 ^ re respectively 
the numerical coefficients of N? and JVf, and whose values change with the indie, al and edge conditions. 

In light of equations (50) through (53), equation (48) could be rewritten as 


EE 

k t 


+ P mn kt + fimnkt 


A U =0 


(54) 


12 


where 


Mmnkt - 


ab 

f?A 66 <x 


0 > mnk£ 


xy 


+ 


-tnnki ( a mnk£ a mnk£ a mnk£ a mnkd + a mnkt ( a mnk£ a mnk£ ~ a mnkt a rnnkl) 


a 22 , . fl 33 _ n 23 „32 

u mnk£ a mnk£ a mnk£ a mnk£ 


P mnkt — 


ab 


r)A 66 a. 


*v 


[^(m, fc)(j 4 iiaa; + A^oty) + C(n, £)(A2iat x + A22&y)] 


(55) 

(56) 


In equation (54), M mnki and P m nkl terms contain material properties which are temperature depen- 
dent. Thus, m the eigenvalue solution process using equation (54), one has to assume a temperature T a 
and use the corrresponding material properties as inputs to calculate the eigenvalue 1/T cr where T cr is 
the buckling temperature. This iteration process must be continued until the assumed temperature T a 
approaches the buckling temperature T cr . 

Thus, in thermal buckling, when the buckling temperature is to be calculated instead of the buckling 
load factors, the eigenvalue solution process requires a temperature iteration process, and therefore, is 

slightly different from that in mechanical buckling for which only a one-step eigenvalue solution process is 
needed. 


6.3 Different Face Sheets Temperatures (Nonzero Thermal Moments) 

When the face sheets temperatures are different (i.e., T x \ T 2 ) the sandwich panel will be subjected 
not only to thermal forces { Aj , Nj , Nj y } but also thermal moments {Mj, Mj, Mj }. The problem will 
t en become a “bending” problem, and no longer an “eigenvalue” problem. One can then calculate the 
panel deflection w in terms of Fourier coefficient A mn . The buckling condition will correspond to the 
condition for which the term in series representation of w, (associated with a particular buckling mode 
shape), becomes unbounded (i.e., A mn — > oc for a given {m,n}). 

As will be seen later, the buckling conditions for the case with thermal moments turned out to be the 
buckling equation for the case without thermal moments. Therefore, only the case with 4S edge condition 
will be analyzed as an example. 


6.3.1 Thermal Moments 


Let the thermal moments {AfJ , , Mj y } be expressed in double Fourier series in accordance with 

the deformation functions given in equations (24) to (26) for the 4S edge condition as 


OO OC 


m= 1 n=~ 1 

OO OO 

<=EE H 

m— 1 n~ 1 

OO OO 


. thtcx . nny 
mn sin - sin 


a 

. mnx . 
l sin sin 




mnx 

mn COS COS 


b 

nny 

~b~ 

nny 


m=l n= 1 


where the Fourier coefficients F mn ,// mn ,5 mn are given by 

rb 


F - — ( a f b \aT ■ mnx ■ nny 

mn ab J 0 J 0 Mx sm a Sm —tedv 


(57) 

(58) 

(59) 

(60) 
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, 4 f a f b ., T . mirx . rnry 

H ™ = al l 

„ 4 f a [ b . - T twtx rnry 

S "" = at l L < T" 5 ' 


(61) 

(62) 


In light of the deformation functions (eqs. (24) to (26)) and the thermal moment expressions 
(eqs. (57) to (59)), the energy equations (8) (setting u = v = 0) and (9) may be integrated to yield 
the forms given in Appendix E. Notice that in equation (E-l), the thermal moment terms turned out to 
be linear functions of {A mn , B mn , Cmru} an d n °t quadratic functions of {A mn , B mn ,C mn }- 


6.3.2 Nonhomogeneous Equations 

After the application of the Rayleigh-Ritz method according to equation (36), one obtains three non- 
homogeneous simultaneous equations for each indicial set of {m, n} shown in the following 


fl mnmn^mn + 32 ^ X/ ^rnnkl^kl 4” a mnmn^Lnn + a T nnrnn^rnn Jn 


„21 


in 4“ Q mnmn^n “1" fl mnmn^n Kmn 


23 


a 


31 


32 

mnmn ^mn 4" a mnmn ^mn 


Bmn “1“ ^mrimn^-Lrin 3 J r 


(63) 

(64) 

(65) 


where the coefficients a^ nmn are defined by equation (C-l), and the nonhomogeneous terms J mn , A mn , 
and Lmn are defined as 


/m7r\ 2 f mr\ 2 /m7r\/n7r\ 

J mn = Fm n (— J + H mn 25 mn ( G ) ( f, ) 


Kmn = -F, 


ran 


f run 


- H (—\+S 

= ~ i- *Jmn 

\ a J 


ri7r\ o fnn\ 

-)+Sm»( T ) 

/n7r\ 

It; 


(66) 

(67) 

(68) 


Prom equations (64) and (65), B mn and C m n may be solved in terms of {A mn , K m n , Lmn, a mnmn 1 as 

.33 _ T „23 

(69) 


o2?_„_a2 l 


a 2 J a? 3 


n u mnmn u mnmn “mnmn^mnmn a 

Bmn ~ a 22 a5? - a 2 J a 32 


21 _ 

„ _ a mnmn a mnmn u mnmn u mnmn a 

(-'mn — 22 „33 _ fl 23 „32 

^mnmn^mnmn u mnmn u mnmn 


33 

'mn 

32 


mn + 22 „33 _ „23 „32 

mnmn u mnmn ~ u mnmn u mnmn a mnmn u mnmn mnmn mnmn 


^ n 33 _ r 

** mn a rnnmn ^mn^mnmn 


a 2 ? a? 1 


+ 


I n 22 

iJmnU’r 


^mn a mnmn 


„32 

a mn ^ 


mn%nmn 


mn T fl 22 


mnmnumnmn u m nmn u mnmn 


33 


- al 3 _a 32 


Substitution of equations (69) and (70) into equation (63) yields, for each indicial set {m, n}. 


Mn 


+ ^2 ^2 fimnktAkt = Rr 


(70) 


(71) 


n xy k t 

where M mnmn is defined in equation (49) under the condition rn=k, n=£, 7 / =32, and R m nmn is defined by 


_ 1 Q k / J _ a mlmn(^m n a mn mil LmnOm n 

mnmn 32fc zy D* \7T / mn °mnmn a mnmn — a mnmn a mnn 


nmn 
mnmn 


) 


+ 


13 /fe' _32 _ 7 „22 \ 

a mnmn\ iv *"" Q mnmn Lmn , wnmnJ 
22 /jo3 

a mnmn a mnmn u mnmn u mnmn 


(72) 
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Equation (71) forms an infinite number of nonhomogeneous simultaneous equations each of which is 
associated with a set of {m n } values (or mode shape) for the calculation of infinite number of Fourier (or 
j coefficients A mn in the series representation of panel deflection w(x, y ) equation (24). 

6.3.3 Buckling Condition 

The calculated Ritz coefficients A mn have the following form 


^mn — 


~ (73) 

nt th ° ™ me *' ator ( contains {M mnmn , R mnmn , <5 mn fc<}, and the denominator A is the determinant 
r p , COC lcien s ° kf - ° 1 e nonhomogeneous simultaneous equations written out from equation (71) 

r"': a (D 12 is shown in Appcndlx D for either <*" < d - , » - 

I he mathematical meaning of the buckling state in light of equation (73) is that the Ritz coefficient 

.."I TT 0 ! U Un v (l f ’ Pand deflecUon - or A ^ 0). Namely, when the buckling state is 

c icd, the term in the series (eq. (24)) which corresponds to the particular deformation mode shape 
becomes the most important term. * 

From the above analysis, one sees that the buckling conditions for the cases with and without the 
thermal moments are identical (i.e., A = 0) under the classical small deformation theory. Because of this 
ing, similar bending analyses for nonzero thermal moments for other edge conditions were not carried 


7. NUMERICAL EXAMPLES 

7.1 Physical Properties of Sandwich Panels 


The sandwich panel is assumed to be fabricated with titanium face sheets and titanium 
core, having the following geometrical and material properties. 


honeycomb 


Geometry: 

h = 1.2 in. 
t s = 0.032 in. 

a = Q o = 24 in. (for varying b ) , or ab = a 2 0 (for constant panel area) 

Material properties: 

Face sheets 



O 

O 

1000 °F 

E x = E y , lb/in 2 

16 x 10 6 

10.5 x 10 6 

G xy , lb/in 2 

6.2 x 10 6 

4.7 x 10 6 

Vxy = Vyx 

0.31 

0.31 

a x = ay, in/in-°F 

4.85 x 10~ 6 

5.6 x 10 -6 

a x y in/in-°F 

0 

0 

p Ti lb/in 3 

0.16 

0.16 
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Honeycomb core (properties at 600 °F) 


Ecx 

— 

2.7778 x 10 4 

lb/in 2 

E C y 


2.7778 x 10 4 

lb/in 2 

Ecz 

= 

2.7778 x 10 5 

lb/in 2 

Gcxy 

= 

0.00613 

lb/in 2 

Gcyz 

= 

0.81967 x 10 5 

lb/in 2 

G C zx 

— 

1.81 x 10 5 

lb/in 2 

v Cxy 

= 

0.658 x 10- 2 


VCyz 

= 

0.643 x 10" 6 


VCxz 

= 

0.643 x 10~ 6 


Otx = oc y 

— 

5.37 x 10 -6 

in/in-°F 

a X y 

= 

0 

in/in-°F 

Phc 

= 

3.674 x 10 -3 

lb/in 3 


7.2 Buckling Interaction Curves 

In generating the data for plotting the buckling interaction curves, k y was set to zero. For a given aspect 
ratio different values of k x were assigned; then the corresponding eigenvalues l/k xy were calculated 
from equation (48). The buckling interaction curves plotted in k x -k xy space for different aspect ratios 
(a = constant) are shown in figure 6. For the case of the square (£ = 1) panel, the additional set of 
buckling interaction curves shown in broken curves is for the case when the effect of transverse s ear is 
neglected. For the square panel, the buckling interaction curves for 4S, 2S2C, and 4C cases are continuous 
curves of symmetric buckling. However, for the 2C2S case, the buckling interaction curve is a composite 
curve, partly for symmetric buckling and partly for antisymmetric buckling. Notice that the effect of the 
transverse shear is quite large for the sandwich panel. Without the consideration of the transverse shear 
effect, the sandwich panel buckling strength could be overpredicted considerably. 1 he 4S case has the lowest 
buckling strength. Clamping the two opposite edges (i.e., from 4S case to 2C2S and 2S2C cases), coul 
enhance the buckling strength considerably. By additional clamping of the other two opposite edges (i.e 
from 2C2S and 2S2C cases to 4C case), further improvement of the buckling strength could be achieved. 
However, the improvement is not as large as that for the previous case (i.e., from 4S case to ' ‘ U1 
2S2C cases). With or without the consideration of the transverse shear effect, the improvement of buckling 
strength through edge clampings is larger in pure compression than in pure shear. 

At aspect ratio £ = 0.7 (a = constant), only the buckling interaction curve for the 2S2C case is 
continuous, and for symmetric buckling only. For the remaining three cases (i.e., 4S, 2C2S, and 4C cases), 
the buckling interaction curves are composite curves, each consisting of two segments: one segment for 
symmetric buckling, the other segment for antisymmetric buckling. 

At higher aspect ratios (i.e., £ = 2,3,4), most buckling interaction curves are composite curves. Some 
interaction curves contain more than two segments. A maximum of four segments could be found for the 
cases of 4C and 2S2C at £ = 4. At high aspect ratio (i.e., 5 = 4) the buckling strengths of the 4S and 
2C2S cases are relatively close, and also those for the 4C and 2S2C cases. 


7.3 Buckling Curves for Pure Compression and Pure Shear 

Figures 7 and 8, respectively, show k x (pure compression) and k xy (pure shear) plotted as functions of 
aspect ratio £. In changing £, the length of the panel (a) was kept unchanged (i.e., a = 24 in.). Notice 
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Omt the rate of reductions of k x and k xy values with the increase of | are the most severe in the region 

fficr'Vt’ n an ,? ® rad ^ a11 ^ d ! crei f e “a increases. At higher aspect ratios (| > 2), the buckling curves of fc. 

4Cea.<L h XV ^ f ° r l ff G 2C2S ^ 2S2C C8SeS ’ respective1 ^ converge toward those for the 4S and 
©cause e edge effects from the two sides (i.e, two vertical edges in figs. 7 and 8) of the panel 
diminish as ~ increases. 

Figures 7 and 8 are the conventional plots of buckling curves. They may not serve as ideal design 
curves for aerospace structural panels because when J is changed (under a = constant), the panel weight 
th ’ t PanG ar f a ° . 1S f° chan S ed accordingly. In the aerospace structural designs, the main objective is 

? ameiyi for a siven panei weight ’ the ***" * ** * 
ptimum buckhng strengths (or stiffnesses). For this reason, modified buckling load factors k x and k xy 

(k y - 0) were recalculated as functions of £ under the condition ab = a 2 0 = constant (instead of a = 

constant). Figures 9 and 10, respectively, show the alternative plots of k x vs £ and k xv vs £ when the panel 
area a* was kept unchanged. In practical applications, the panel has *to be%upporL by frame 

(cross section assumed constant), therefore, the edge frame weight (or edge frame length, (a + b)/2a 0 ) 
was also plotted in figures 9 and 10 as a function of £. The square panel (£ = 1) has the minimum ed E e 
frame weight; however, it does not have the optimum buckling strengths in compression and shear. The 

— W bUCklmg strengths ( fig - 9 ) reached minimum at £ = 1.6, 1.4, 2.2, 1.0, respectively, for the 4S, 
C, 2C2S, and 2S2C cases. The lowest shear buckling strengths (fig. 10) occur at £ =09 09 12 07 
respectivdy, for the 4S, 4C, 2C2S, and 2S2C cases. Figures 9 and 10 serve as design cJ’s trsllLing 
the desired sandwich panel geometry (i.e, £ value). To boost the panel buckling strengths in compression 
and shear, some weight penalty due to edge frame is inevitable. The desirable hifh-stif Lss-LZ^Z 
panel shapes will be slightly slender (£ < 1). " 

7.4 Buckling Interaction Surfaces 

ed Jrnl!iMB°'Tn^r e 'fi di ^ mensional buckling surfaces which are plotted in {k x , k x „, i) space for different 
j \ , e 1 ® ure > domains of symmetric and antisymmetric buckling (lowest buckling; 

odes) are also shown. Figure 1 1 gives better visualization of the buckling behavior of the sandwich panel 

tan the traditional buckling plots shown in figures 6 to 8. For slender rectangular panels (i e 4 < 1 ) 
antisymmetric bucklings occur mostly in the compression-dominated regions. For wider panels (i e 4 >’ 

t'hLkl L?K ym ur etriC ’7 kli " SS If 11 ' P,aCe in the ®hea.r-dom ina ted regions. In the neighborhood of 4 - 1 

Onl!lr ~ .7, Tu m t k ' ™ - *• » - >• *> fOT “>e 4S, 4C, and 2820^.' 

y oMhe 2C2 case s the lowest buckling mode in the compression-dominated region antisymmetric 

( 771 > n » “g. 4). Such buckling behavior also occurs in the flat rectangular plates of £ ss 1. 

7.5 Thermal Buckling Curves 

For most of the practical materials, the coupling coefficient of thermal expansion a, „ is zero Therefore 
generating the data for thermal buckling curves, the following conditions were imposed. 

= <* y , a xy = 0 (i.e, 7Vj y = Mj y = 0) (74) 

Ml = = 0 (75) 

The above conditions will induce m-plane biaxial thermal compression without shear and bending 
igure 12 shows the critical buckling temperature (T cr ) plotted as a function of |, with the panel length a 
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being kept constant. Those thermal (biaxial compression) buckling curves somewhat resemble the uniaxial 
compressive buckling (mechanical buckling) curves shown in figure 7. For buckling temperatures higher 
than 1000 °F, the face sheet material property data at 1000 °F were used as inputs to equation (54) for 

calculations because of the lack of material property data at high temperatures. For the honeycomb 
core, the only available material property data at 600 °F had to be used as inputs for T cr calculations. It 
was found that T cr was relatively insensitive to the material property change with temperature. 

For the present particular panel (i.e., dimensions chosen) the thermal buckling temperatures T cr exceed 
the titanium melting point (3074 °F) at low J, and gradually decrease with the increase of At high 
aspect ratios, for the 4S and 2C2S cases level off at about 1000 °F (below superplastic temperature, 
1650 °F), and for 4C and 2S2C cases, at temperatures slightly below the melting point. 

Figure 13 shows the alternative plots of T cr as a function of £ for constant-area panels (i.e., ab° = 
constant) the lowest buckling temperatures for 4S, 4C, 2C2S, and 2S2C cases are, respectively, 1297 F, 
3702 °F (above melting point), 2194 °F (above superplastic temperature), and 2205 °F (above superplastic 

temperature), and which occur respectively at | = 1.0, 0.975, 1.8, and 0.5. 

For the present sandwich panel, the actual thermal buckling will take place only for the 4S case in the 
region 1.5 < J < 1.8. Outside this region for the 4S case and for all the whole range of 7l for the other 
three edge conditions, no “real” thermal buckling could take place because the sandwich panel will first 
undergo superplastic creep or melting. 

8. CONCLUDING REMARKS 

Through the use of the Rayleigh- Ritz method of minimizing the total potential energy of a structural 
system, the combined load (mechanical or thermal) buckling equations were established for orthotropic rect- 
angular sandwich panels under four types of edge conditions. Twodimensional buckling interaction curves, 
and three-dimensional buckling interaction surfaces were constructed for high-temperature honeycomb-core 
sandwich panels supported under four different edge conditions. The buckling interaction surfaces provide 
easy visualization of the variation of the panel buckling strengths, and the domains of buckling modes 
(symmetric and antisymmetric) with the edge condition. Furthermore, the buckling temperature curves 
for the sandwich panels were presented. 

The effect of transverse shear on the buckling strength is quite large for sandwich panels, and by 
neglecting the transverse shear effect, the buckling strengths could be overpredicted considerably. With the 
inclusion of the transverse shear effect, the buckling load factors became panel-size-dependent in addition 
to panel-aspect-ratio dependent. Clamping the edges could enhance the buckling strength greatly more in 
compression than in shear. Thermal buckling conditions for the cases with and without thermal moments 
were found to be identical for the small deformation theory. 
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Figure 1. Honeycomb-core sandwich panel. 
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Figure 2. Combined compressive and shear loadings of a rectangular sandwich panel. 
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Figure 3. Deformation of a sandwich panel in the x-z plane. 
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Figure 4. Four types of edge conditions. 




Figure 5. Forces and moments acting on the differential element of 


a sandwich panel. 
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Figure 6 . Combined load buckling interaction plots for honeycomb-core sandwich panel under different 










( c ) a ~ 2 > a = 3. | = 4 respectively. 


Figure 6. Concluded. 






Figure 7. Comparison of compressive buckling strengths of honeycomb-core sandwich panels under different 
edge conditions, a — constant. 
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Figure 9. Comparison of compressive buckling strengths of honeycomtvcore sandwich panels under different 
edge conditions. Constant panel areas. 
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Figure 10. Comparison of shear buckling strengths of honeycomb-core sandwich panels under different 
edge conditions. Constant panel areas. 





Figure 11. Buckling interaction surfaces for honeycomb-core sandwich panels under different edge condi- 
tions. a = constant. 
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(c) Two sides clamped, two sides simply supported (2C2S). 
Figure 11. Continued. 





Figure 12. Thermal buckling temperatures for honeycomb-core sandwich panels under different edge 
conditions, a = constant. 
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APPENDIX A 
INTEGRAL RELATIONS 

In the integrations of the potential energies V, and V 2 , the following integral relations were used. 
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For m = k \ 1 : 
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For m = k or m ^ k: 
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For m ± k — odd or even: 
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APPENDIX B 
ENERGY EQUATIONS 

After substitutions of the trial deformation functions (eqs. (24) through (35)) into the energy equations 
for V\ (eq. (11)) and V 2 (eq. (9)) (signs of forcing functions reversed according to eq. (12) or (13)), and 
carrying out double integrations using the special integral relations given in Appendix A, the energy 
components AV\ and AV 2 respectively for V’i and V 2 have the following expressions for different edge 
conditions under particular indicial conditions. 

Case 1. Four edges simply supported (4S edge condition). 
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Case 2. Four edges clamped (4C edge condition). 
(1) m = k ^ 1, n = t = 1 
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(4) m = Jfc k 1, n = l % 1 
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(7) to - fc = 2, w = l = 1 


** - See 


- ^ (1 + fc) 4 + Dqg ^ (1+A;) 2 A m i - - B k \ 

- 2 (D 12 + D 66 ) (|) (1 + *) 2 ] (^A ml - f* ml ) (^A kl - C fcl ) 


4D22 (?) +d “GD (I+t); 


£i4 ml -CmlHj4bl-CM 


H" fc) Bm\Bk\ — DQyCmiCkl 


^=|EE{[^(;) 2 a + *) 2 + ^(f) 2 ]^ tl } 

( 8 ) m - fc = 2 ,n=^^l 

A* 2 . = £ E e : e {- 1 k (z) 2 (> + *) 4 + n- m 2 0 + ‘) ! (> + ^ 


m n - 2 k 


^ f ^ ^mn Brrinj ( ^ ^4/cn 


- (£>12 + £> 66 ) ^ — j (1 + k) 2 ( 1 + n 2 )J - Bmn'j (^ A kn ~ C kn j 

— 2 ^ 22 [(1 + n 2 ) 2 + 4 n 2 ] + £>66 ^ (1 + k) 2 (l + n 2 ) (—A mn — C n , 


^ ^ Akn Ckn 


■(1 + k) 2 B mn B kn -£83L(l + n 2 ) C mn C kn 


AK 2=§EEE{^ ^(^(l + fcf + W^^l+n 2 ) A mn A k 


41 



(9) m — fc = 2 , n — 1 = 2 
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Case 3. Two sides clamped, two ends simply supported (2C2S edge condition). 
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Case 4. Two sides simply supported, two ends clamped (2S2C edge condition). 
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APPENDIX C 

COEFFICIENTS OF CHARACTERISTIC EQUATIONS 

The coefficients a l ^ nkt of the simultaneous characteristic equations (37) through (39) are defined in the 
following for different edge conditions under different indicial conditions. 

Case 1. Four edges simply supported (4S edge condition). 
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Case 2. Four edges clamped (4C edge condition). 
(1) m=n=k=t=l 
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(5) m = fc = l, n — 1 = 2 
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a mnmt = a mnmt = ^ |^22 (^j + if + (£*12 + 2-Dg6) ^ ^ m2 )(l + ^) 2 | 

^ — 2 |^ U (a) (l + m2 ) 2 +4 m2 + £>66 (1 + m 2 )(l + if + Dq x {\ + m 2 )| 


(C-7) 


n 22 

u mnmt 


„23 _ „32 _ _ Im j. n \ /nX/7T 

® J mnmt O'mnmt 0 {*-'12 i ^66/ 

^ V a 


-^ ) (1 + m 2 )(l + if 


a 33 = _ I 

U 'mnvnt 2 


D22 {^) (l+^) 4 + Ai6 Q (1 +m 2 )(l + £) 2 + D Qy {\ + tf 


(C-8) 
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°mlfcl — ~ 


(7) m - k = 2,n = l= \ 

[ D n (-) (1 + k) 4 4- 2(£>i2 + 2D 66 ) j ^ (l + £) 2 +4 D 2 2 

+ D * ( 5 ) a [!** ( 5 ) a f 1 + *> 2 + *v (f ; 

“mlfcl = a mltl = \Du (1 + k) 4 + (Di 2 + 2 D m ) (1 + k) : 

“Zki = a mifci = 4D 22 (£)* + ( D 12 + 2Z? 66 ) (1 + k ) 2 

i^ 11 (a) ( 1 + *0* + -^66 (1 + k) 2 + ^Dq x (1 + ky 


n 22 - - 


a mlfcl — a mUl — ~ (^12 + ^6e) (1 + k ) 





„33 - 

a mlfcl — 


4D 22 (^) + A>e(^) (1 + k ) 2 + D Qy 
(8) m - A: = 2, n = ^ ^ 1 


(C-9) 


®mnfcn 


— — 2 |^ n (a) (1 + k) 4 + 2(D 12 + 2D 66 ) ^ ^ (1 + A:) 2 (l + n 2 ) 

+ D 22 (^) 4 [(1 + n 2 ) 2 + 4n 2 ] } + ^ (^) 2 [fcx (^) 2 (1 + k) 2 + k y (^) 2 (1 + n 2 ) 


n 12 _ n 2l _! 

Tnnfcfi mnfcn 2 


„13 - n 31 - ± 

TTitikn mtikn 2 


Dn © (1 + fc)4 + (z?12 + 2Dfi6) © ©* ( ! + fc ) 2 (! + 


n 


n 22 - - I 

u mnfcn 2 


^22 
2 



1 + n 2 ) 2 + 4n 2 + (Z?i 2 + 2/^66) ( — 



(1 -+- A:) 2 (l + u 2 ) 


£>u © ( 1 + /c ) 4 + D 66 © (l + fc) 2 (l+n 2 ) + £> Qx (l +A:) : 


1 


°mnkn — a mnkn = ~ 12 + £>66) ! (1 + fc) 2 (l + n 2 ) 



„33 - _ I 

mnfcn 2 


j^M© [(1 +n 2 ) 2 +4n 2 ] + D 66 © (1 + fc) 2 (l + n 2 ) + D Qy (l + n 2 ) j (C-10) 
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(9) m-k = 2 ,n-£ = 2 


,11 


a mnkt — J 


£>n fo) (1 + k) 4 + 2 (Di 2 + 2 Dqg) f-') ' (l+k) 2 (l +£) 2 


+ £>22 


(t) 4 c^ 



„1 2 _ _21 __1 
a mnkl Q'mnkt ^ 

„13 _ -31 _ 1 


^ ll (o) + ^) 4 + (£>12 + 2£> 66 ) ^ j (1 + fc) 2 (l + £) 2 

D22 (f ) (1 + l) A + (£>12 + 2 D 66 ) f-) 2 (f) (1 + fc) 2 (l + £) 2 


22 J- 

®mn« = ^ 


Du (£) (1 + fc)4 + Dm (f ) 2 (1 + fc)2(1 +*) 2 + d q *<- 1 + k ) 5 


a mnM — iw = 7(^12 +£>66) ( — ) (— ) (1 + fc) 2 (l + if 


£>22(1) (l+^) 4 + £>66 (-) (1 + fc) 2 (l + if + D Qy (\ + £f 


n 33 _ I 

TTlTlki. ^ 


Case 3. Two sides clamped, two ends simply supported (2C2S edge condition). 
(1) m = h, n = £ = 1 

.4 / \ 2 / _ \ 2 




_3 Al (v) +8(D 12 +2 A»)(™) (£) +16£> 22 (£) 


A2 


-D*{- 

v a 


21 


( m7T\ 



V 1 ml a mlml ~ 


„13 = „ 3 * - _ 

a mlml a mlml 


Skx \~) + 4 Mt 

3£>u (v) + 4 ( D 12 + 2D 66 ) ( 


run \ ( 7 r 


\ a 



16fl22 ^) 3 +4(fl , 2+2flM) (^) 2 Q- 


«mlml = 3 £*1 1 ) + 4£>66 f 2 + 3 Dq 


23 _ 32 

a m 1ml ^mlml 


2 33 

x m 1ml 


= 4(D, s + D 65 ) (™) (I) 

= 16£> 22 + AD m {^pj 2 + AD Qx 


(C-ll) 


(C-12) 
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(2) m = k, n = l ^ 1 


Omnmn = 2D n + 4(D i2 + 2D 66 ) (jj (1 + n 2 ) + 2 D 22 [(1 + n 2 ) 2 + 4 n- 

-^)>(v) 2+ ‘«(f) 2 <H 

-’“"(”) ,+s t D “ +sn - ) (") (?)*(* + "*>] 

= - |2I?22 (j) [(1 + n 2 ) 2 + 4n 2 ] + 2(D X2 + 2 D 66 ) (1 + n 2 ) | 


a 12 = a 21 

^mnmn ^mnmn 


a 13 = a 31 

^mnmn u mnmn 


.22 


= 2£>ii(™J +2A»(£) (l + n 2 ) + 2D Qx 
= 2(D i2 + D 66 ) (j) (1 + n 2 ) 


a 23 = a 32 

u 'mnmn ^mnmn 


11 

a mnm£ 


°mnmn = 2^22 [(1 + n 2 ) 2 + 4n 2 ] + 2D 66 (1 + n 2 ) + 2D Qy (l + n 2 ) 

( 3 ) m = k, n — 1 = 2 

£>n (^) 4 + 2{Dn + 2D ™ } ( j ) 2 (1 + 1)2 + ° 22 ( i ) 4 (1 + ^ 

^)>(™) 2+ *»a) 2 H 

a mnmi = <£w = D U + Pl2 + 2D 66 ) (^) (1 + *) 2 

“mnm< = “mnm< = £>22 (j) (1 + ^) 4 + (£>12 + 2 D 66 ) (^) (* + *)* 

£>11 ( V) 2 + ^ (f ) 2 (1 + £)2 + ° QX 


n 22 = 


“mnml = “mnm< = ■ (£>12 + £> 6 g) ( 5 ) ( 1 + 

2 2 

£> 22 ^^ (1 + £) 4 + £>66 (1 +£) 2 + Dq v (\ + ty 


„33 

a mnmi 


(C-13) 


(C-14) 
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Case 4. Two sides simply supported, two ends clamped (4S4C edge condition). 
(1) m = k = 1, n = t 


a lnin = 16D U + 8 (D l2 + 2 Z? 66 ) (£) + 3D 22 


V \b 


«i»i» = o?i ln = - 16D U (^) +4(D 12 + 2D 66 )^) ^ 
®lnln = ®lnln = “ 3D 22 ^ ^ + 4(D 12 + 2 D 66 ) (^~ 


«Sm = I 6 O 11 (J) S + 4Z* (^)%4D 
«?nln = ^In = 4(D 12 + £> 66 ) Q ) 


a lnln — 3Z) 22 


+ 4£> 66 - + 3Dr 


(2) m = k ^l,n = l 


a — ^ 

^mnmn — a mnmn t 


a mnmn - 2 ^11 [(1 + 7B 2 ) 2 + 4m 2 + 4(D 12 + 2£> 66 ) j (^) (1 + m 2 ) + 2L>22 j 

- 2D 'G) 2 Kf) 2 < 1 + ”* 2 > + *»(T)l 

a mnmn = °mnmn = - |2£>11 j (1 + m 2 ) 2 + 4m 2 J + 2 (£>i 2 + 2£>66) ("y') 0+ m2 )j 

a mnmn ~ °mnmn t 5 — |^D 22 ^ + 2(Di 2 + 2Z?66)(1 + Tn?) ^ 

fl mnmn = 2D u Q [(1 + m 2 ) 2 + 4m 2 ] + 2£> 6 6 (1 +m 2 ) + 2D Qy (l +m 2 ) 

“mnmn = «“ „ mn = 2(D 12 + D 66 ) (l) (H) (1 + m *) 


mnmn u 'mnmn 


-2 £> 22 ^ +2£) 6 6^j (1 + m 2 ) + 2Dq v 


(C-16) 
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(3) m — k — n = £ 

= - Dm (;)V + *=)“ + 2(»i* + 2D 66 ) (^) 2 (^) 2 (1 + k) 2 + D-n (™) 

+ D ‘(^) k *{a) (1 + * )2+ * : »( t ) ] 

aJU = <£!„*„ = Dn (~) 0 + k)> + (Du + 2Dee) (l) ( y)’ (1 + *) 2 
“mnta = “m„kn = D 22 (y) + (D, 2 + 2D 66 ) Q (t) ^ + ^ 

“mnkn=- (~) 0 + *) 4 + £>M ("£•) (' + k ) 2 + D Qi (\ + k) 2 
= <4 2 „*„ = - (Dn + Dee) (^) ( ™) 0 + *) 2 

4?„l„ = ~ D 22 (^) 2 + Dee(^) 2 (l+t) 2 +D 0! , <°- 17 > 
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APPENDIX D 
BUCKLING EQUATIONS 

The buckling equations (eigenvalue solution equations) written out from equation (48) up to order 12 
(i.e., 12 x 12 matrices) for the case m±n = even (symmetric buckling) and m±n = odd (antisymmetric 
buckling) for different edge conditions are given below. 

Case 1. Four edges simply supported (4S edge condition) (ref. 9). 
m±n = even (symmetric buckling) (4S) 


Akt — 1 > 


A13 

^22 

M 1 

Ai* 

A 24 

^33 

A 42 

Asi 

^35 

A44 

^53 

m = l,7i=l 

fc X y 

0 

4 

9 

0 

0 

8 

45 

0 

8 

45 

0 

0 

16 

225 

0 

m=l,n=3 


^1313 

K X y 

4 

5 

0 

0 

8 

7 

0 

8 

25 

0 

0 

16 

35 

0 

m=2,n=2 



M222 

fciy 

4 

“5 

20 

“63 

0 

36 

25 

0 

20 

"63 

4 

7 

0 

4 

7 

7 \ 

ii 

6 




Ml 31 

k Z y 

0 

8 

“25 

0 

8 

7 

0 

0 

16 

35 

0 

m = l,n=5 





Ms 15 

k X y 

40 

"27 

0 

8 

"63 

0 

0 

16 

"27 

0 

m=2,n=4 






A*2424 

k X y 

72 

“35 

0 

8 

63 

8 

3 

0 

120 

147 

m— 3,n— 3 



Symmetry 



M333 

^xy 

72 

35 

0 

0 

144 

49 

0 

rn=4,n=2 








M242 

kxy 

40 

27 

120 

“147 

0 

8 

3 

171=5,71=1 









MSl 151 

0 

16 

“27 

0 

m=3,n=5 










k X y 

80 

21 

0 

771=4,71=4 











M 444 

kxy 

80 

21 

m=5,n=3 












AfttH 


k 

K xy 
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m ± n = odd (antisymmetric buckling) (4S) (ref. 9) 


A-kt — * 

A\2 

A21 

Ma 

^23 

A 32 

A 41 

A IQ 

A 25 

^34 

A43 

^52 

A61 


Maui 

Kxy 

4 

0 

4 

0 

_ _8_ 

0 

20 

0 

8 

25 

0 

4 

_ ^ 

m=l,n=2 

“9 

5 

45 

63 


%50 

m=2,n=l 


A^2121 

fciy 

8 

”45 

0 

4 

5 

0 

4 

“35 

0 

8 

25 

0 

20 

63 

0 

771— 1 ,77=4 



Afl414 

kxy 

8 

"7 

0 

16 

”225 

0 

40 

27 

0 

16 

"35 

0 

8 

“ 175 

m -2,n=3 




Af 2.72.7 

kxy 

36 

”25 

0 

4 

”9 

0 

72 

35 

0 

4 

“ 7 

0 

m=3,n=2 





Af.72.72 

^xy 

8 

“7 

0 

4 
" 7 

0 

72 

35 

0 

4 

9 

77i=4, n= 1 






Af4141 

K X y 

8 

“ 175 

0 

16 

“35 

0 

40 

27 

0 




0 




Af 1616 

kzy 

20 

0 

8 

0 

36 

to 

II 

Ti 

E 



oyinineiry 



”11 

45 

1225 

m=2,n=5 








Af2525 

^xy 

8 

“3 

0 

100 

“441 

0 

771 = 3,77 = 4 









Af3434 

k X y 

144 
“ 49 

0 

8 

"45 

m=4, n=3 










Af 4343 
fc X y 

8 

”3 

0 

m=5,n=2 











Afs252 

fc x y 

20 
" 11 

777 = 6 , 77 = 1 












Affilfil 

fcxy 


(D-2) 
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Case 2. Four edges clamped (4C edge condition). 
m =L n = even (symmetric buckling) (4C) 


Akt —> 



^22 

a 3 i 

^15 

-424 

-433 

A 42 

Asi 

A 35 

-444 

^53 

m=l,n = l 

^1111 

k X y 

Afl H3 
fC X y 

4 

225 

Muai 

kxy 

0 

8 

” 1575 

Afi IM 

k X y 

8 

“ 1575 

0 

0 

16 

11025 

0 

m=l,n=3 


A^1313 

fcxy 

44 

” 1575 

Ml Mi 

k, X y 

Miiii 

k X y 

184 

4725 

Mini 

k X y 

88 

11025 

0 

Af 1335 

k X y 

368 

33075 

0 

771=2,77=2 



M2221 

fc X y 

44 

1575 

4 

525 

M 2224 

fcxy 

484 

11025 

A^2242 

fc X y 

4 

525 

44 

“3675 

M221± 

kxy 

44 

”3675 

m=3j7i=l 




-^3131 

k x y 

0 

88 

11025 

^3133 

fc X y 

184 

4725 

^3151 

fc X y 

0 

368 

33075 

M 3153 

kxy 

771 — 1,71=5 





fcxy 

104 

“2079 

M 1533 

kxy 

8 

3675 

0 

Mi535. 

kxy 

208 

14553 

0 

771=2,71=4 






Af242i 

fc X y 

2024 

“33075 

Af2442 

kxy 

8 

“ 3675 

104 

1323 

M 2441 

^iy 

184 

11025 

771=3,71 = 3 



Symmetry 



Af 3.133 

k x y 

2024 

”33075 

Afaasi 

k xy 

M333& 

fc X y 

8464 

99225 

M2 3S3 
kxy 

771=4,71 = 2 








^4242 

k X y 

104 

“2079 

184 

11025 

Mi211 

fc X y 

104 

1323 

771=5,71=1 









M 5151 

k X y 

0 

208 

14553 

Af-S153 

kxy 

771 = 3,71 = 5 










AT3535 

k X y 

4784 

43659 

A^3553 

kxy 

771 = 4,71 = 4 











MlUl 

fc X y 

4784 

43659 

771 = 5,71 = 3 












M 5353 

kxy 


(D-3) 
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m±n = odd (antisymmetric buckling) (4C) 


Akt -» 

A 12 

•^21 

A14 

•<423 

^32 

- An 

Al6 

^25 

A 34 

A 43 

^52 

A 61 

m=l,n=2 

Afl212 

*xy 

4 

A*1214 

K X y 

44 

1575 

Afi232 

fc X y 

8 

0 

4 

Afl234 

*xy 

88 

0 

4 


225 

1575 

“525 

" 11025 

4725 

m=2,n=l 


A^2121 

fc X y 

8 

■^2123 

i& 

“4725 

44 

1575 

Afel41 

iT 

4 

0 

88 

M 2143 

4 

0 



1575 

4725 

“ 11025 

A: X t/ 

"525 



A^1414 

^xy 

M 1432 

Muis 

104 

Afl434 

36& 

0 

8 

m= 1 ,n~ 4 



^xy 

“11025 

^xy 

2079 

k X y 

33075 

" 33075 





A*2323 

484 

•^2341 

*xy 

172 

A^2325 

2024 

M2343 

44 

0 

m— 2,n— 3 




fcxy 

" 11025 

17325 

fcxy 

33075 

fcxy 

3675 






-^3232 

184 

0 

44 

^3234 

2024 

Mi252 

172 

m=3,n— 2 





h X y 

“4725 

3675 

k X y 

33075 

kxy 

17325 







M *141 

8 

0 

368 

Ml43 

K X y 

104 

Mlfil 

m=4,n=l 






k X y 

" 33075 

33075 

2079 

fc X y 




Symmetry 




236 

^1634 

344 

0 

4 

771 = 1 ,77=6 





foxy 

"3861 

k X y 

" 121275 

“ 99225 









^2525 

104 

Mto43 

4 

0 

m—2,n=5 








k X y 

1323 

k xy 

" 1225 










Mi 434 

8464 

Ml452 

344 

II 

y? 

II 

E 









kxy 

" 99225 

k X y 

" 121275 

m=4, n=3 










Mama 

104 

^4361 










k X y 

” 1323 

k X y 

m=5,n=2 











^5252 

k X y 

236 

"3861 

771=6, n=l 












MSl ifiJL 

k X y 


(D-4) 
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Case 3. Two sides clamped, two ends simply supported (2C2S edge condition). 
m ± n = even (symmetric buckling) (2C2S) 


Aki — > 

An 

An 

M 2 

Mi 

^15 

Ma 

Mz 

M 2 

Mi 

^35 

A 44 

M3 

m~ 1 ,n=l 

Mini 

fcxy 

Min i 

k X y 

4 

45 

0 

0 

8 

“315 

0 

8 

225 

0 

0 

16 

" 1575 

0 

771=1,71=3 


Mim 

kxy 

44 

"315 

0 

M 1315 

kxy 

184 

945 

0 

88 

" 1575 

0 

0 

368 

4725 

0 

771=2,71=2 



M2222 

k X y 

4 

“25 

4 

105 

A^2224 

kxy 

44 

175 

0 

4 

"63 

12 
" 175 

0 

44 

441 

m=3,n= 1 




M 3131 

k xy 

0 

8 

175 

M 3 133 
kxy 

8 

35 

0 

0 

16 

245 

0 

m=l,7i=5 





MlSlS 

kxy 

520 

"2079 

0 

8 

525 

0 

0 

208 

"2079 

0 

771 = 2,71=4 






A*2424 

kxy 

184 

"525 

0 

8 

441 

104 

231 

0 

184 
" 1323 

m =3,71=3 



Symmetry 



M 3323 

kxy 

88 

"245 

0 

Maaas. 

kxy 

368 

735 

0 

771=4,71=2 








M 4242 

K X y 

8 

"27 

24 

245 

M 4244 

kxy 

88 

189 

771 = 5,71= 1 









Mmsi 

kxy 

0 

16 

189 

M 5153 

k X y 

771=3,71 = 5 










Masas 

k 

*-zy 

1040 
" 1617 

0 

771 = 4,71=4 











M4444 

kxy 

368 

"567 

771=5,71=3 












M 5353 

k 

A 'xy 


(D-5) 
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m±n — odd (antisymmetric buckling) (2C2S) 



A 12 

Mi 

A\4 


A 32 

A 4 I 

Al 6 

A 25 

A 34 

A 43 

^52 

Agi 

m=l,n=2 

Aft 212 

K X y 

4 

25 

Aft214 

K X y 

44 

315 

0 

8 

225 

0 

4 

105 

0 

88 

1575 

0 

4 

175 

m=2,n=l 


Aftl21 

*xy 

8 

315 

Aft 123 

^xy 

4 

25 

0 

4 

925 

0 

8 

“175 

0 

4 

65 

0 

m=l,n=4 



Mini 

^iy 

184 

945 

0 

16 

1575 

Ml416 

k X y 

520 

2079 

0 

368 

“4725 

0 

8 

1225 

m— 2,n=3 




A/2323 

k 

^xy 

44 
~ 175 

0 

172 

3465 

K X y 

184 

525 

0 

44 

“441 

0 

m— 3,n=2 





A^3232 

k X y 

8 

35 

0 

12 

175 

Af3234 

^xy 

88 

245 

0 

4 

45 

m=4,n=l 






Aft Ml 

*xy 

8 

4725 

0 

16 

245 

Aft 143 

K X y 

8 

27 

0 

m=l ,n=6 



Symmetry 



k X y 

1180 

3861 

0 

344 

17325 

0 

4 

3675 

m=2,n=5 








A42525 

fcxy 

104 

231 

0 

4 

147 

0 

m=3,n— 4 









Aft434 

fcxy 

368 

735 

0 

8 

315 

m=4,n=3 










Mimi 

kxy 

88 

189 

0 

m=5,n— 2 











A/5252 

k X y 

4 

” 11 

m=6,n=l 












MfiM 
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Case 4. Two sides simply supported, two ends clamped (2S2C edge condition), 
m ± n = even (symmetric buckling) (2S2C) 


Akt —> 

A n 

-4l3 

-422 

^31 

-4i5 

-424 

-433 

-442 

-4 5 i 

-435 

A 44 

-4 5 3 

m~ l,n— l 

M ini 
fcxy 

0 

4 

45 

7^7 1131 
k X y 

0 

8 

225 

0 

8 

315 

0 

0 

16 

" 1575 

0 

m = l,n=3 


Afl313 

4 

“25 

0 

0 

8 

35 

M 1333 

k X y 

8 

175 

0 

0 

16 

245 

0 

m=2,n=2 



A/2222 

fcry 

44 

“315 

4 

“63 

0 

44 

175 

A/2242 

kxy 

4 

105 

44 

441 

0 

12 
“ 175 

m— 3,n— 1 




Mini 

k X y 

0 

88 

1575 

0 

184 

945 

M 3151 

^iy 

0 

368 

4725 

0 

m=l,n=5 





Afl515 

K xy 

8 

“27 

0 

8 

441 

0 

A/ 1535 
kxy 

16 

189 

0 

m= 2,n=4 






A /2424 

k X y 

88 

“245 

0 

8 

525 

88 

189 

A /2444 

kxy 

24 

245 

m=3,n=3 



Symmetry 



A/arm 

fc X y 

184 

“525 

0 

0 

368 

735 

A/,1353 

kxy 

771—4,71—2 








A*4242 

k X y 

520 

“2079 

184 

1323 

0 

104 

231 

771 = 5,71=1 









Msisi 

fcxy 

0 

208 

2079 

0 

771 = 3,71=5 










A*3535 

kxy 

368 

“567 

0 
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m ± n = odd (antisymmetric buckling) (2S2C) 
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APPENDIX E 

ENERGY EQUATIONS WITH NONZERO THERMAL MOMENTS 

Using a similar process mentioned in Appendix B, the energy equations with nonzero thermal moments 
will be presented in the following for the 4S edge condition as an example. 
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where the last term of equation (E-2) is nonzero only under the indicial conditions: m ± k = odd and 
n ± 1 = odd. 
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